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Abstract 

The mixed phase of layered superconductors with no magnetic screening is studied through a 
partial duality analysis of the corresponding frustrated XY model in the presence of weak random 
point pins. Isolated layers exhibit a defective vortex lattice at low temperature that is phase 
coherent. Sufficiently weak Josephson coupling between adjacent layers results in an entangled 
vortex solid that exhibits weak superconductivity across layers. The corresponding vortex liquid 
state shows an inverted specific heat anomaly that we propose accounts for that seen in YBCO. 
A three-dimensional vortex lattice with dislocations occurs at stronger coupling. This crossover 
sheds light on the apparent discrepancy concerning the observation of a vortex-glass phase in recent 
Monte Carlo simulations of the same XY model. 
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High-temperature superconductors are layered and extremely type-IlQ. This fact led to 
the suggestion early on that a state with phase-coherent yet decoupled layers is possible Q]. 
It was later demonstrated, however, that any amount of Josephson coupling between layers 
results in a macroscopic Josephson effect across layers at low temperature [jj. In the presence 
of external magnetic field perpendicular to the layers, other workers made the analogous 



proposal that sufficiently weak coup 
two-dimensional (2D) vortex lattices 



ine, could lead to a decoupled stack of phase-coherent 



i 41 J. Monte Carlo simulations!™ and a partial duality 

an 

analysis [7] [8] of the corresponding frustrated XY model demonstrate, however, that a highly 
entangled vortex lattice state with relatively small or no phase coherence across layers does 
not exist in practice. 

The elusive decoupled vortex lattice state may exist at low temperature in the presence 
of random point pins, however jg||lol|. In this paper, we show that this is indeed the case 
through a partial duality analysis of the corresponding layered XY model with uniform 
frustration^]. We show first that a defective vortex lattice state can exist in isolated 
layers. It exhibits macroscopic phase coherence in the zero-temperature limit, despite the 
presence of unbound dislocations that are assumed to be quenched in by the random pins 12]. 



We next turn on Josephson coupling between adjacent layers and 



ind that weak supercon- 



1]| in the zero-temperature 



ductivity exists across layers at sufficiently high layer anisotropy 
limit. After assuming a continous 2D ordering transition for each layer in isolation, we then 
find that an inverted specific heat jump can occur inside of the vortex liquid state at weak 
coupling. This prediction compares favorably with the recent observations of such a peak in 
the high-temperature superconductor YBCojl^]. 

2D. Consider a stack of isolated superconducting layers in a perpendicular external mag- 
netic field. In the absence of Josephson coupling as well as of magnetic screening, the XY 
model over the square lattice with uniform frustration provides a qualitatively correct de- 
scription of the mixed phase in each layer. The corresponding Boltzmann distribution is set 
by the sum of energy functionals 



e ( xUi) = - E E h cos ( A ^ - 4J (!) 



,(2) 

\ b ) — - 

for the superfluid kinetic energy of each layer I written in terms of the superconducting 
phase 0(r, /). Here A Ai 0(r, I) = 0(r + a/t, I) — 0(r, /) and A = (0, 2nfx/a) make up the local 
supercurrent, where / denotes the concentration of vortices over the square lattice, with 



lattice constant a. The local phase rigidity J^if, I) is assumed to be constant over most of the 
nearest-neighbor links (r, r + afi) in layer I, with the exception of those links in the vicinity of 
the pinning sites that are located at random. After taking the Villain approximation, which 
is generally valid at low temperature Q| , a series of standard manipulations then lead to a 
Coulomb gas ensemble with pins that describes the vortex degrees of freedom on the dual 
square lattice[8]. The ensemble for each layer / is weighted by the Boltzmann distribution 
set by the energy functional 

EUl) = (27T) 2 £ 6Q J G< 2 > 5Q' + £ V p \Q\ 2 , (2) 

(R,R') R 

written in terms of the integer vorticity field Q(R,l) over the sites R of the dual lattice in 
that layer, and of the fluctuation SQ = Q — f '. A logarithmic interaction, = — V~ 2 , 
exists between the vortices, with a strength J equal to the gaussian phase rigidity Last, 
V P (R, I) is the resulting pinning potential [8|]. 

We shall next assume that the array of random pins in each layer, V P (R, I), quenches in 
unbound dislocations into the triangular vortex lattice at zero temperature [12j. To check for 
superconductivity in such a defective 2D vortex lattice, we now compute the macroscopic 
)hase rigidity, which is given by one over the dielectric constant of the 2D Coulomb gas (j2J) 
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pf D) /J = 1 - Iim(27r/77 SW ) (SQ^Q^/^a 2 ^ . (3) 



Here SQ^ = — (Qt) is the fluctuation in the Fourier transform of the vorticity in layer 
I'- Qk = Q{R"> l) e%k R - Also, ?7 SW = kBT/2iiJ is the spin-wave component of the phase- 
correlation exponent, and A/ii denotes the number of points in the square-lattice grid. Now 
suppose that a given vortex is displaced by Su with respect to its location at zero tempera- 
ture. Conservation of vorticity dictates that its fluctuation is given by —V • 5u. Substitution 
into Eq. (jnj) then yields the result 

pf^l Jo = 1 - WJrhw) (4) 
for the phase rigidity in terms of the vortex component to the phase-correlation exponent, 

R 

The latter monitors fluctuations of the center of mass of the vortex lattice [16(. Above, iV vx 
denotes the number of vortices, while a vx = a/ f 1 ^ 2 . 



To proceed further, we now express the displacement field as a superposition of pure wave 

n 

and defect components of the triangular vortex lattice|l6j]: 5u = 5u wv + Su^i- Notice by Eq. 
(JHJ) that phase coherence is insensitive to the wave contribution if rigid translations are 
excluded, since J2^^wv = in such case. The former is achieved through bulk pinning^ |, 
and the latter then follows under periodic boundary conditions. Consider now a single 
unbound dislocation with Burgers vector b that slides along it glide plane yj] a distance 5Rdf 
with respect to its location at zero temperature. The relative displacement field, 5Hdf, then 
corresponds to that of a pure dislocation pair of extent SR d f that is oriented along its glide 
plane. After following steps similar to those taken in ref. [16] for the pristine case, it can be 
shown that Eq. (|SJ) yields a fluctuation of the center of mass 



V ' vx = n d{ (\6R d{ \ 2 )(b/2a vx y\nR /a Ai (6) 

for the vortex solid, where n d { denotes the density of unbound dislocations, where adf is 
the core diameter of a dislocation, and where Rq is an infrared cut-off. Above, the overbar 
denotes a bulk average. Observe now that both rj sw and r]' vx vanish linearly with temperature. 
By Eq. (0J), we conclude that the defective 2D vortex lattice shows a positive phase rigidity 



in the zero-temperature limit at sufficiently dilute concentrations o 



unbound dislocations, 
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of the 2D Coulomb 



n di — > 0. The above is borne out by direct Monte Carlo simulations 
gas ensemble (J2J). 

The previous positive result for macroscopic phase coherence [Eq. (0}] in the zero- 
temperature limit can be confirmed by calculation of generalized phase auto-correlation 
functions within an isolated layer: Ci[q] = (exp[i J2r ' 0(^O])o- Following a similar 
calculation in the pristine case[16|, application of the Villain approximation [see Eq. (J2J) 
and ref. [14(] yields the form Ci[q] = \Ci[q\\exp[i Y^rQ^^oi^ 0] f° r these autocorrelations, 
where (f>o(r,l) represent the zero-temperature configurations of isolated layers. In the low- 
temperature regime, phase correlations are then found to decay algebraicly like 

\d[q)\ = 9o + ■ exp[ V2D £ g(l)ln(r 12 /r ) q(2)} (7) 

(1,2) 

at the asymptotic limit, r 12 — > oo, with a net correlation exponent approximately equal to 
t] 2 d — Vsw + v'vx f° r small vortex components, rj' vx ?7 SW . Here, g = p@ D ^ / J is the ratio of 
the 2D stiffness with its value at zero temperature, J, while n + counts half the number of 
probes in q(r). Also, r denotes the natural ultraviolet scale. 



3D. We shall now add a weak Josephson coupling energy — J 2 cos(A z — A z ) to all of the 
verticle links in between adjacent layers of the three-dimensional (3D) XY model. Here, 
J z = J/7' 2 is the perpendicular coupling constant, with anisotropy parameter 7' > 1. 
The layered XY model can be effectively analyzed in the selective high-temperature limit, 
fcgT 3> J z , through a partial duality transformation. This leads to a dilute Coulomb gas 
(CG) ensemble that describes the nature of the Josephson coupling in terms of dual charges 
that live on the vertical links. Phase correlations across layers can then be computed from 
the quotient 

;exp[i$>(O0(O]) = Z CG [p}/Z CG [0} (8) 



of partition functions for a layered CG ensemble |1 11]: 

Zcg[ P }= E Vo^Mah] ■ e-^ r ^ % (9) 
{*».«} 

where the dual charge, n z (f, I), is an integer field that lives on links between adjacent layers 
I and I + 1 located at 2D points r. The ensemble is weighted by a product of phase auto- 
correlation functions for isolated layers I probed at the dual charge that accumulates onto 
that layer: 

Qi(r) = P(r, I) + n z (f, I - 1) - n z (f, I). (10) 

It is also weighted by a bare fugacity yo that is raised to the power iV[n z ] equal to the total 
number of dual charges, n z = ±1. The fugacity is given by y = J z /2ksT in the selective 
high-temperature regime, J z ksT, reached at large model anisotropy. 

In the absence of Josephson coupling, random point pins lead to zero-temperature phase 
configurations, (po(f, I), that are completely uncorrelated across layers. At zero parallel field, 
Eqs. (JBJ and Q therefore yield the expressions |3] 



(cos0 M+1 ) - y o £a(0, 1) ■ C7 +1 (0, 1) (11) 
1 

and 



l(e i0i ' !+1 )l 2 = 3/o EE C ^ ^(O, 2) • a+i(0, 2)Q +1 (0, 1) (12) 
1 2 

for the inter-layer "cosine" and the inter-layer phase correlation, to lowest order in the fu- 
gacity. The overbar represents a bulk (disorder) average, while 0^ + i(f) = (f)(r,l + 1) — 
0(r, I) — A z (f, I) is the gauge-invariant phase difference across adjacent layers. Macroscopic 
phase coherence shown by each layer in isolation (jJJ) is lost at a transition temperature [l7| 



Tj 2D ) > 0. Only short -range phase correlations on the scale of ^2D exist at higher 
temperature following C/(l,2) = ^"^/^e^^e - ^ ^. By analogy with 2D melting 
physics|3|Q|) the presence of quenched-in unbound dislocations also implies that only 
short-range phase correlations exist inside of each layer in isolation, on average, at zero tem- 
perature. Specifically, we have exp[i0^ +1 (l)] • exp[— z0| ( |^ 1 (2)] = e~ 2ri2 ^ 2D asymptotically, 
where 0^ +1 (r) = 4>o(r, I + 1) — <po{r, I) — A z (r, I) is the quenched interlayer phase difference, 
and where 1 2 d represents a zero-temperature disorder scale set by n d{ . Substitution into ex- 



pression (O then yields the result | (e^w+i) p ~ [g%( J/k B T)(l 2D £ 2D / Al)] 2 for the inter-layer 
phase correlation inside of the critical regime, ^> ho, where A = j'a is the Josephson 
penetration length. This approximate result reaches unity at a cross-over field 

/ 7 'x 2 ~ g 2 (J/k B T)(l 2D £ 2D /al), (13) 

in units of the naive decoupling scale $o/Aq, that separates 2D from 3D vortex-liquid 
behavior Q. Substitution into expression (jllj) for the inter-layer "cosine", on the other 
hand, yields a non-divergent result 

<cos0 M+1 > ~ g 2 (J/k B T)[(l£ + ^y'/Ao} 2 (14) 



that is va' 



be shown 



id in the decoupled vortex liquid that exists at fields much larger than f"f'd- It can 



1 2(| that the next-leading-order term for the inter-layer "cosine" is negative, 
that it diverges just like the leading order term for the inter-layer correlation ()12|). and that 
it becomes comparable to its own leading order term precisely at fields below the 2D-3D 
cross-over scale, Eq. (fT3*j) . Last, Eq. (pPfj) implies an anomalous inter-layer contribution to 
the specific heat per volume equal to 

6($ - 2[1 + ^ 2D /l 2D )]-\d\n^/dT)ej, (15) 



where ej = (cos0^ + i) ■ J/A^d is the Josephson energy density, and where d denotes the 
spacing in between adjacent layers. It also notably shows an inverted specific heat jump that 
is followed by a tailoff at a temperature T p such that ^ 2 d{T p ) ~ 1 2 d if 62D diverges faster 
than (T — T^ 2D ^) _1 . This approximate result is again valid at high anisotropy, 7' > 7^, 
which yields the bound Z 2 d < Q^iksT / J) 1 ^ 2 Aq on the 2D disorder scale by Eq. (JT3J). 

The previous analysis clearly demonstrates that a selective high-temperature expansion 
in powers of the fugacity y necessarily breaks down in the ordered phase, T < T^ 2D \ where 



<^2D is infinite. At this stage it becomes useful to re-express the layered CG enesemble © by 
replacing Ci[q] with its magnitude (J7J), and by compensating this change with the additional 
replacement of A z {f,l) with — (jffl, -, (r). A Hubbard- Stratonovich transformation of the CG 



partition function reveals 
(LD) model with an energy functional that is given bv jlll 

1 



2l| that it is equivalent to a renormalized Lawrence- Doniach 



^LD = P f D) 



J i 



2 (W,) 2 -A W M+1 1 (16) 



where #n+i = 4>fi + i + Oi+i — Oi- A standard thermodynamic analysis[7|] then yields that the 



IT] 



strength of the local Josephson coupling is given by (cos0^ + i) = yo + go(cos9ij + i). It can 



also be shown 



111 ] that phase coherence exists across a macroscopic number of layers, with 



a corresponding phase rigidity equal to pjr / J z = <7o(cos 



In order to compute (cos^+i) at low temperature, we must first determine the configu- 
ration that optimizes -Eld- Eq. (fTHjl implies that it satisfies the field equation 

- V 2 ^ (0) + A - 2 [sin 9[% - sin flgj = 0. (17) 

In the weak-coupling limit, Aq — > oo, we therefore have that 9\ \ if) is constant inside of each 



layer. The fact that e^We - *^ 00 ) = then implies that cos#^ + i = at zero temperature 
in the weak-coupling limit. Indeed, the LD "cosine" can be calculat ed p erturbatively, where 
one finds that cos6^/ + i ~ (/ 2 D/Ao) 2 ln(Ao//2Z)) at zero temperaturejl(J| j^. In the opposite 
limit of weak disorder, 1 2 d — * oo, Eq. (|T7|) yields that sin^j^ l = sin 9ii +1 , on the other hand. 
This then implies that cos^^i = 1 in the weak disorder limit. The bulk average cosfyj+i 
at zero temperature must therefore pass between zero and unity at Ao ~ ho- This condition 
defines a decoupling cross-over field /7^(0) ~ {ho/^y) 2 in units of $o/Aq, at which point 
the reversible magnetization shows a broad diamagnetic peak. By the discussion following 
Eq. ()16|) , we conclude that random point pins result in a vortex glass at sufficiently high layer 
anisotropv 3 [lol] . A ^> 1 2 d, that exibits weak superconductivity across layers [l^: pj: -C J z . 

The results of the above duality method are summerized by the phase diagram displayed 
in Fig. n The present theory notably predicts that an inverted specific heat anomaly (jump 
followed by a tailoff) occurs at weak coupling in the vortex liquid when the 2D correlation 
length £ 2 d matches the 2D disorder scale 1 2 d if %2D diverges faster than (T— J^ 21 ^) -1 [see Eq . 
(fT5|)]. Such a feature has in fact been observed within the vortex liquid phase of YBCO 13]. 
The weight of the latter peak is about Ae exp — 6mJ/cm 3 , while the peak shown by Eq. 



([To]) has a weight Aej = ($o/167t 3 A|Aq) A(cos 0z,z+i). Equating these and using values of 
A L = 140 nm and A = 7nm for the penetration depths in YBCO0 yields a 10 % jump in 
the "cosine". 

Last, although recent Monte Carlo simulations of the same XY model studied here do 



m n 2 = s, 



indeed find evidence for a phase-coherent vortex glass at|24j f^' 2 = 16 and at 
another one using f^y 12 = 2 does not. We believe that the zero-temperature cross-over 
shown in Fig. ^ between an entangled vortex glass and a 3D vortex lattice containing 
dislocations is the origin of this discrepancy. 

In conclusion, a duality analysis of the XY model finds that random point pinsj^fiol] 
drive a cross-over transition in the zero-temperature limit between defective vortex lattices 
that show strong versus weak superconductivity across layers as a function of the Josephson 
coupling[ll]. We further propose that the inverted specific heat anomaly observed recently 
inside of the vortex liquid phase of YBCO does not signal a phase transition |13|, but rather 
is due to the thermodynamic resonance found here, Eq. (fl3|) . 

The author thanks O. Bernal and Y. Nonomura for discussions. 

Note added: Recent Monte Carlo simulations of the same XY model studied here also 
find a non-critical specific-heat anomaly in the vortex-liquid phase (see ref. 
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FIG. 1: Shown is the proposed phase diagram assuming weak point pins and a continuous vortex 
glass phase transition for isolated layers. The concentration of in-plane vortices, /, is held fixed, 
and a mean-field temperature dependence, J oc T c q — T, is assumed. Monte Carlo simulations of 
the same XY model studied here find evidence for a second-order transition between the vortex- 
glass and the vortex-liquid phases (see ref. 



22|). When confronted with the first-order decoup ling 
transition that is expected to separate the vortex liquid from the 3D vortex lattice (see refs. |rjj 
and ^ ms i m P nes the existence of a critical endpoint consistent with experiments on YBCO 



(ref. |l3j) and with other numerical simulations (ref. |24l| ) o f the present XY model. A transition 
to a 3D vortex lattice without defects is reported in ref. 25] at fj 12 < 1. 
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